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Abstract. We use a Green's function method to study the temperature- 
dependent average moment and magnetic phase-transition temperature of the striped 
antiferromagnetism of LaFeAsO and other similar compounds as the parents of FeAs- 
based superconductors. We consider the nearest and the next-nearest couplings in the 
FeAs layer and the nearest coupling for inter-layer spin interaction. The dependence of 
the transition temperature Tn and the zero-temperature average spin on the interaction 
constants are investigated. We obtain an analytical expression for Tn and determine 
our temperature-dependent average spin from zero temperature to Tn in terms of 
unified self-consistent equations. For LaFeAsO, we obtain a reasonable estimation 
of the coupling interactions with experimental transition temperature Tn = 138 K. 
Our results also show that a non-zero antifcrromagnctic (AFM) inter-layer coupling is 
essential to the existence of a non-zero Tn and the many-body AFM fluctuations reduce 
substantially the low-temperature magnetic moment per Fe towards the experimental 
value. Our Green's function approach can be used to other FcAs-based parent 
compounds and these results should be useful to understand the physical properties of 
Fe As-based superconductors. 
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1. Introduction 

The discovery of high temperature superconductor LaF^FeAsOi-x by Kamihara 
et al[l] has triggered world-wide researches on all aspects of FeAs-based pnictides 
superconductors and their parent compounds, namely LnFeAsO (Ln=La[l, 2, 3, 4, 5, 
6, 7], Ce[8, 9], Pr[10], Nd[ll, 12], Sm[13, 14, 15], ...) and AFe 2 As 2 (A=Ca[16, 17, 18], 
Sr[19, 20, 21], Ba[22, 23, 21]). LnFeAsO and AFe2As2 own some common characteristics: 
a) they are all of layered structures and have structure transitions from high- 
temperature tetragonal to low-temperature orthorhombic symmetry; b) they all have 
stripe-like antiferromagnetic (AFM) order formed by Fe atoms and the AFM transition 
temperatures Tn's are not higher than the structure transition temperatures Ts's; c) the 
onset of superconductivity competes with the AFM and structure transitions [24, 25]. 
And they differ in the sense that Tn < Ts for LnFeAsO and Tn = T$ for AFe2As2- Both 
of the two series can be made superconducting by doping them with appropriate dopants 
or applying pressures. It should help understand the superconductivity of FeAs-based 
materials to elucidate the corresponding antiferromagnetism of the parent compounds. 

LaFeAsO is the prototype and the representative of the parent compounds of 
FeAs-based superconductors and thus we focus on the striped AFM order of undoped 
LaFeAsO, whose Tn is 138K and Ts 156K[2, 3]. First-principles results confirm that the 
stripe-AFM order is the magnetic ground state[5]. Spin-wave approaches were adopted 
to give the low-temperature excitation spectra[7, 26, 19], and the spin-orbit interaction 
and p-d hybridization are used to understand the observed small magnetic moment 
0.25yU B ~ 0.36yU B per Fe at low temperature [2, 3, 6]. However, it is highly desirable to 
describe the magnetic moment from zero temperature to Tn within a unified theory. 

In this paper, we use a Green's function method to study the temperature- 
dependent average moment and phase-transition temperature of the striped 
antiferromagnetism of LaFeAsO and other similar compounds as the parents of FeAs- 
based superconductors. We consider the nearest and the next-nearest couplings in the 
FeAs layer and only the nearest one for the inter-layer spin interaction. The dependence 
of the transition temperature Tn and the zero-temperature average spin on the four 
interactions are investigated. We obtain an analytical expression for Tn and determine 
our temperature-dependent average spin from zero temperature to Tn in terms of unified 
self-consistent equations. For LaFeAsO, we obtain a reasonable estimation of the 
coupling interactions with experimental phase-transition temperature Tn = 138 K. Our 
results also show that a non-zero antiferromagnetic inter-layer coupling is essential to the 
existence of a non-zero Tn and the many-body AFM fluctuations reduce substantially 
the low-temperature magnetic moment per Fe towards the experimental value. More 
detailed results will be presented in the following. 

The remaining part of this paper is organized as follows. In next section, we shall 
give our spin model, the Green's function derivation and our main analytical results. In 
section 3, we shall present our numerical results and make corresponding discussions. 
And our conclusion is given in section 4. 
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2. Effective Model, Green's function derivation, and main analytical results 

To deal with the striped AFM configuration of LaFeAsO, we consider the Fe lattice 
of the original orthorhombic LaFeAsO structure and divide it into two sublattices in 
each of which the Fe spins align parallel but antiparallel between the two sublattices 
(figure 1). Hence we consider an anisotropic Heisenberg Hamiltonian 

^ = E J W)Si-S j + J 2 ^S i -S j , (1) 

<ij> 

in which Si denotes the quantum spin operator at the lattice position i, (i,j) means 
nearest-neighbour (NN) spin pairs, and ((i,j)) means next-nearest-neighbour (NNN) 
spin pairs in a-b plane (we only consider NNN pairs in a-b plane because the inter-layer 
interactions are very weak). NN interaction can be three values: J\ a which is the 
spin interaction between parallel NN spins in a-b plane, Jib between antiparallel NN 
spins in a-b plane and J c between inter-layer NN spins. J 2 is interaction between NNN 
spins in a-b plane. Four different J's make H anisotropic. To differentiate spin operators 
in SLI and SL2, we use Su and S 2 j to represent them respectively. For spins in SL2, we 
make transformations: S^j = ~&2ji = ^2j> = ^2j an d then have 

Su ■ S 2 j = -(S'liS'y + •S'li'S'y) + s^s^ 

I (2) 

_ / Q+ QI+ I Q— qi— \ QZ C/Z 

~ 2^ li 2 J ~ t ~ °li°2jJ °li° 2 j • 

Inserting (2) into (1) we can get the Hamiltonian expressed by Sii and S' 2 j (simple but 
too long to give out here). 




Figure 1. The magnetic cell with volume a x 2b x 2c of the orthorhombic Fe spin 
lattice. The Fe lattice consists of two spin sublattices: SLI (black) and SL2 (red or 
gray), a and b are the two base vectors in the FcAs layer, and c is perpendicular to 
both of a and b. 

Accordingly, we use Green's function method [27, 28] to solve the model (1). In this 
scheme one uses a double-time Green's function ((A(t); B(t'))) (A and B represent two 
arbitrary quantum operators) which satisfies the following equation of motion: 

ih±({A(t);B(t'))) = 

5(t-t')([A(t),B(t)\) + ({[A(t),H];B(t'))). (3) 
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This approach proves successful for various Heisenberg spin models [27, 28, 29, 30, 
31, 32]. In [30], there is a detailed but somewhat long derivation about the Green's 
function method for Heisenberg spin model. Accordingly, we sum up the derivation 
in [30] and get the point that the process of using Green's functions to solve the 
average spin of Heisenberg model can be simplified into three steps: a) construct Green's 
functions S^)) and their equations of motion via (3) and then use Tyablikov cutoff 
approximation (4) to decouple the equations of motion[28]; 

{(S?S^B))^{S Z ){(S^B)) (4) 

b) use spectrum theorem to express the correlation function (S~S + ) in term of the 
average spin ^-component (S z ) and then get the Q((S Z )) function 

mS z }) = liS-S^/iS*) ; (5) 

c) use the Callen Expression (6) [30] to evaluate (S z ) self-consistently 

~ (s - m + ®) 2S+1 + (s + 1 + m 2S+1 ff , 

^ ' (1 + $)2S+1 _ $2S+1 • V°) 

According to the three steps given above, for our spin system we construct double- 
time spin Green's functions between spin operators at two positions i and j in SL1 

Gg 1) (t,f) = «^[(*);^i(0» (7) 

and Green's functions between spin operators at two positions i' in SL2 and j' in SL1 

Gg?(*,O = <0£<O;^(O» • (8) 

Equation (7) can be expressed as Fourier expansion 



G£ 1) Me- iw(t - t ' )/B dw, 



(9) 



2nh I 

because the Hamiltonian (1) is time independent, and (8) is similar. Assuming that 
each spin has the same average (S*) for SL1 and (S'2) for SL2 and because of the 
AFM symmetry plus the transformation (2), we have (S*) = (S^) = (S z ). Then using 
the Fourier transformation as shown in (9), making Tyablikov cutoff approximation to 
decouple the equations of motion and another Fourier transformation from lattice sites 
real space to k space, we can have 



(11) 



•AftPi(k) 



UJ 



and 



in which 



2{S Z )± 

^k 1)j lbP 2 ( k ) +#k 21> J lbPi( k ) 



J lb g^p 2 (k)=0 



UJ 



2(S Z ) 



(10) 



(11) 
2k 



(21) 



r 



uj e 



-ikr 



-ikr 



(12) 
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p 1 (k) = —p(l — cosk-a,) + 1 + 2q + r 
p 2 (k) = (1 + 2q cosk-a) cosk-6 + r cosk-c 



(13) 



and 



P 



(14) 



Jla J2 Jc 

— , q = — , r = — . 

■Jib Jib Jib 

We should point out that: a) the wave vector k we used here is based on the whole 
lattice sites (SL1+SL2), so r in summations of (12) runs over all sites in the whole 
lattice; b) for homogeneous system Gfl] r is only a function of relative position r and 
independent of i (as a result G^ = G {11) which is used below); Gf^ r is analogous. 
From (10) and (11) we derive 



(11) 



(S z ) XPi + yfp^Pi Px-y/pFp.. 



\Ip\-p\ 



to 



E(k) 



to 



E(k) 



where E(k) is the spin excitation spectrum defined by 



E(k) = 2J lb (S z ) x /pl(k)-pi(k). 
And from <j£ we get G~ 



G s i) M4£ 3relk ' N) 



(15) 



(16) 



(17) 



keBZ 



in which N is the total number of spins in SL1 and SL2, and BZ denotes the first 
Brillouin zone (there are N k-points in BZ). Using spectrum theorem and letting j = i 
we get the correlation function (S^ Sf) as follow 



(SiSt) 



1 

N 



Im[G (11 V + zO+)] 



E 

keBZ 



Pi 



duo 



P\~P\ 



coth 



;is) 



where (3 = l/(k B T), T is temperature, and k B is the Boltzmann constant. Then using 
(5) we get the <3> function: 



keBZ V ^1 



Pi 



coth 



f3E(k) 



Pt 



(19) 



Now, the average spin ^-component (S z ) can be obtained easily by self-consistently 
solving (19) and (6). A special case is that when the temperature T = 0, coth — > 1 
and we have 



T=0 



1 

2iV 



E 

keBZ 



Pi 



(20) 



At this time, $0 is no longer dependent on (S z ) and the zero-temperature average spin 
^-component (S Z ) Q can be obtained directly by insert (20) into (6). 
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While temperature approaches to T N , (S z ) approaches to zero and further E(k) — > 
and $ — > oo. Expanding (19) and (6), we derive 



<S Z >°M/1-|-, (21) 



where T N is defined by 

m 2j lb s(s + r 



3Tk B 

andr = iE k [P 1 /(P?-P?)]- 

3. Numerical Results and discussions 



(22) 




(b) (c) (d) 

*r v *f- 



Figure 2. (a) Structure of FeAs layer and scheme for the exchange interactions 
mediated by Fe-As-Fe paths: Fel-Asl-Fe2 and Fel-As2-Fe2 for J la , Fe2-As2-Fe3 
and Fe2-As3-Fe3 for Ju,, and Fel-As2-Fe3 for J^- Different AFM configurations are 
shown: (b) stripe-AFM along a direction, (c) stripe-AFM along b direction, and (d) 
checkerboard AFM. Energies of the three AFM configurations are given in (23). 

About the coupling interactions J ia , Ju, J2 and J c in FeAs-based pnictides, there 
is no consensus on their magnitudes and many authors only consider two or three of 
them [4, 5, 33, 23, 19, 16, 34]. Yildirim's first-principles results show that J\ ~ ^[5]- 
However, we prefer the opinion that J\ a , Ju, and J 2 originate from AFM superexchange 
through As atoms[35, 36, 37]. From the viewpoint of the structure of FeAs layers, 
both Ji a and Ju, are mediated by two Fe-As-Fe paths (figure 2(a)) and there should 
be Ji a 7^ Ju but Ji a ~ Ju due to the small structure change from tetragonal to 
orthorhombic symmetry. J 2 is mediated by only one Fe-As-Fe path (figure 2(a)) and 
there should be 2J 2 ~ Jia- From the viewpoint of classical favorable energy to form 
stripe-like AFM patterns along a direction as shown in figure 2(b) (see figure 2(b)-(d) 
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Figure 3. The reduced AFM transition temperatures (T^/[JibS(S + 1)]) as functions 
of q for r — 0.001 and p = 0.90, 0.92, • • ■ , 1.0 (from top to bottom). The inset shows 
an r dependence of Tn for p = 0.98 and q = 0.5. 



and (23)), 



E(b) 


= 4J la 5 2 


- u lb s 2 - 


iJ 2 S 2 


E (c) 


= -4Ji a 5 2 


+ AJ lb S 2 - 


iJ 2 S 2 


< E (d) 


= -4Ji a S 2 


- 4J lb S 2 + 


SJ 2 S 2 




< E (c) => 


J\a < Jib 






< E (d) => 


2J 2 > J\ a 





there should be J la < J lb and 2 J 2 > Ji a , which in fact are just the conditions that fulfill 
p 2 — p 2 ^ to make E(k) in (16) meaningful. As for J c , it's a very weak long-range 
AFM interaction with a nowadays unclear origin other than superexchange. Therefore, 
in terms of p, q, r in (14) we confine the coupling interactions as follow: 

1 ^ p > 1 - 6 p 

p/2<q<: p/2 + 5 q , (24) 
< r ^ 5 r 

where 8 P , 5 q , and 5 r are the confine parameters for p, q, and r respectively and they 
fulfill < S p , S q , <5 r < 1. Lhs of (24) represents the necessary conditions to form striped 
AFM ordering along a direction, and rhs the conditions to limit p, q and r within small 
regions. 

From figure 3 we can see that Tn increases as q and r increase but decreases as 
p increases. It's easy to understand. NNN spins and inter-layer NN spins all align 
antiparallel, so bigger AFM coupling interactions will lower the system's energy, stabilize 
the AFM configuration and hence enhance Tn; on the contrary, spins along a direction 
align parallel but have AFM interactions, therefore bigger p will increase the system's 
energy, destabilize the AFM configuration and hence decrease Tsf. It's notable that 
while q — > p/2 or r — > 0, then Tn — * 0, that is to say, both the existence of an AFM 
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inter-layer interaction J c and the condition 2 J 2 > J\ a are essential to form striped AFM 
ordering. In fact, when r = this system becomes two-dimensional, and the result of 
Tn = in two dimensions is analogous to the Mermin- Wagner theorem for isotropic 
interactions [38]. We also see that the critical condition p — 1 doesn't lead to Tn — ■> 0, 
which manifests J\ a = Jib isn't a fatal factor to kill Tn but only a critical value to 
separate the two cases shown in figure 2(b) and 2(c). 

0.8 
0.7 
0.6 



0.4 
0.3 
0.2 

0.45 0.5 0.55 0.6 

q 

Figure 4. The zero-temperature average spins (S z ) as functions of q for r = 0.001 
and different p values: 0.90, 0.92, • • • , 1.0. The lower limits of (S z ) are shown by the 
dotted line. The inset shows (S z ) ~ r curves for (p,q) = (0.98,0.50) (solid) and 
(0.98,0.49) (dotted) respectively. 

Magnetic moment per Fe of LaFeAsO at low temperature is reported as 0.36yUB[2] or 
0.25/ie[3] experimentally, both of which are very small compared with the first-principles 
values 2.2 ~ 2.4/iB/Fe[39, 40]. Accordingly, we choose S = 1 in our model (assuming 
the Lande ^-factor equals 2). Figure 4 shows the zero-temperature average spin (S z ) . 
Similar to Tn, (S z ) is also an increasing function of both q and r but a decreasing 
function of p. However, (S Z ) Q — > only when both r — > and q — > p/2. When only 
r — > or q — > p/2 is met, (5 Z ) approaches to a minimum but not zero, while at the 
same time T N — ► 0. 

There are four J's in our model. What are their values? Then let us have a look at 
what can they be under condition Tn = 138K. Figure 5 shows the regions available for p, 
q and r under conditions (24) with 5 P = 0.05, 5 q = 0.05 and S r = 0.001 in orange (gray) 
colour for J^, = 40meV, 50meV, and 60meV respectively. The smaller J^, the smaller 
the parameter region available. Hence, for given 8 P , S q and S r , there is a lower limit for 
Jib to fulfill a given T N . This lower limit, written as </™ n , is given in table 1, which 
shows that the smaller each of S p , S q and S r is, the bigger J™ n is. However, although no 
upper limit for J u is given out, Jib can not be infinitely great. In fact, first-principles 
results show that Jib ~ 50meV[4, 33]. 

Here we didn't use the experimental data for low-temperature magnetic moment 
per Fe atom, which amounts to {S z ) = 0.13 ~ 0.18, to determine the J's, because 




q=p/2 case r 

j i i 




Figure 5. The parameter regions (orange or gray) satisfying (24) with Tn = 138K 
for Jib = 40meV (a), Ju, = 50mcV (b), and Ju, = 60mcV (c), where we use 
S. p = 0.05, S q = 0.05, and S r = 0.001. The star shows the value (p, q, r)=(0.955, 
0.526, 0.0009) in (a), (0.98, 0.53, 0.0003) in (b), or (0.99, 0.53, 0.0001) in (c). 



it's too small. If {S z ) = 0.18 is met, there have to be 8 r < 10" 6 even if q = p/2 
is taken to minimize {S z ) (see inset of figure 4). Such tiny r definitely cannot fall 
within the orange (gray) area in figure 5 with a reasonable Jib to fulfill Tn = 138K. 
That is to say, although the many-body AFM fluctuations substantially reduce the low- 
temperature magnetic moment per Fe, yet they cannot be in full charge of the very 
small low-temperature magnetic moment which indeed is also ascribed to spin orbit and 
p-d hybridization etc[6]. 



Table 1. 


J^ in for 


given 5 P , 


5 q and S r . 


5 q 


<),. 


s p 


JjfVmeV 


0.01 


0.0001 


0.01 


104.4 






0.1 


89.8 




0.001 


0.01 


87.1 






0.1 


72.6 


0.05 


0.0001 


0.01 


50.9 






0.1 


44.6 




0.001 


0.01 


43.6 






0.1 


37.2 



We choose three sets of J's from the regions available in figure 5 for estimation 




Figure 6. The average spins (S z ) as functions of the temperature T for Jif,=40meV, 
50mcV, and 60meV. The corresponding (p,q,r) parameters are (0.955, 0.526, 0.0009), 
(0.98, 0.53, 0.0003), and (0.99, 0.53, 0.0001), respectively. 



(the three stars): (p, q, r)=(0.955, 0.526, 0.0009), (0.98, 0.53, 0.0003), and (0.99, 0.53, 
0.0001) for J 1{ ,=40meV, 50meV, and 60meV respectively. In terms of J u , J la , J 2 and 
J c , we can refer to them as Ju, = 50 ± lOmeV, J\ a = 49 ± lOmeV, J 2 = 26 ± 5meV, and 
J c = 0.020 ± 0.015meV. The average spin (S z ) vs temperature T curves with the three 
sets of parameters given above are shown in figure 6. 

We take LaFeAsO for example here, however our calculations are not restricted 
to LaFeAsO, because nearly all LnFeAsO has similar even the same transition 
temperatures (see table 2). It seems that doesn't vary much with different Ln. 



Table 2. AFM and structure transition temperatures Tn and Ts of LnFeAsO. 



Ln 


La 


Ce 


Pr 


Nd 


Sm 


T N (K) 


138 


140 


127 


141 


140 


r s (K) 


156 


155 


153 






Ref 


[3] 


[8] 


[10] 


[11] 


[41] 



This is different from v4Fe2As2, whose Tn varies with A: Tn=172.5K[18] for v4=Ca; 
T N =198K[20], 205K[42] or 220K[43] for A=Sr and T N =140K[22] or 135K[44] for A=Ba. 
Although there are some differences of structure between LnFeAsO and AFe2As2, we 
believe that our model works well for both of them, because they both have layered 
structures with stripe-like AFM order formed by Fe atoms. Indeed, this can also be 
extended to Fe(SeTe)[45, 46], which also have nearly the same properties and can be 
superconducting under certain conditions. 
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4. Conclusion 

In summary, we use a Green's function method to study the striped AFM order formed 
by Fe atoms in LnFeAsO which is the representative of the parent compounds of recently 
discovered Fe-based superconductors. We take LaFeAsO for example to analyze the 
AFM transition temperature Tn and zero-temperature average spin (S z ) , and show 
that both T N and (S z ) are increasing functions of Ju, J 2 and J c but decreasing 
functions of J\ a . By using Tn = 138K, we make a reasonable estimation of the 
coupling interactions, Ju, = 50 ± lOmeV, J ia = 49 ± lOmeV, J 2 = 26 ± 5meV, and 
J c = 0.020 ± 0.015meV. Average spin (S z ) is determined in the same way from zero 
temperature to T N and T N is expressed analytically. Our results also show that a non- 
zero AFM inter-layer coupling J c is essential to the existence of a non-zero Tn and 
that the AFM fluctuations substantially reduce the low-temperature magnetic moment 
towards the small experimental value. Although our results cannot determine the 
relations between structure and AFM transitions, we believe that the AFM transition 
is likely caused by the structure transition because most of experimental results show 
Tm ^ Ts for FeAs-based pnictides. Our Green's function approach of the striped AFM 
properties can be used to other FeAs-based parent compounds. 
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